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In  this  paper  we  first  show  that  the  linear  models  of  proton  exchange  membrane  (polymer  electrolyte 
membrane,  PEM)  and  solid  oxide  (SO)  fuel  cells,  commonly  used  in  power  and  energy  literature,  are  not 
controllable.  The  source  of  uncontrollability  is  the  equation  for  pressure  of  the  water  vapor  that  is  only 
affected  by  the  fuel  cell  current,  which  in  fact  is  a  disturbance  in  this  system  and  cannot  be  controlled 
by  the  given  model  inputs:  inlet  molar  flow  rates  of  hydrogen  and  oxygen.  Being  uncontrollable  these 
models  are  not  good  candidates  for  studying  control  of  dynamic  processes  in  PEM  and  SO  fuel  cells. 
However,  due  to  their  simplicity,  they  can  be  used  in  hybrid  configurations  with  other  energy  producing 
devices  such  as  photovoltaic  (solar)  cells,  wind  turbine,  micro  gas  turbine,  battery  (ultra  capacitor)  to 
demonstrate  some  other  phenomena,  but  not  for  control  purposes  unless  the  hybrid  models  formed  in 
such  hybrid  configurations  are  controllable.  Testing  controllability  of  such  hybrid  models  is  mandatory. 
Secondly,  we  introduce  some  algebraic  constraints  that  follow  from  the  model  dynamics  and  the  Nernst 
open-loop  fuel  cell  voltage  formula.  These  constraints  must  be  satisfied  in  simulation  of  considered  fuel 
cell  modes,  for  example,  via  MATLAB/Simulink  or  any  other  computer  software  package. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  controllability  and  observability  concepts  are  the  system 
state  space  concepts.  They  have  been  known  to  control  engineers 
for  more  than  fifty  years  since  the  initial  work  of  Kalman  [1  ].  Slowly 
these  concepts  are  becoming  known  and  used  in  other  engineer¬ 
ing  and  scientific  disciplines,  especially  when  the  so-called  Kalman 
system  canonical  decomposition  was  derived  in  [2,3].  The  Kalman 
canonical  decomposition  states  that  only  the  system  modes  that 
are  both  controllable  and  observably  appear  in  the  system  trans¬ 
fer  function  and  those  either  uncontrollable  or  unobservable  cancel 
out  from  the  transfer  function  (system  input/output  description). 
This  result  has  established  the  fact  that  the  state  space  system 
description  (via  system  eigenvalues)  is  more  general  than  the  sys¬ 
tem  description  via  transfer  function  (via  system  poles)  since  the 
set  of  system  eigenvalues  is  broader  than  the  set  of  system  poles 
(all  the  poles  are  the  eigenvalues,  but  not  all  the  eigenvalues  are 
the  system  poles). 

In  the  first  part  of  the  paper  we  show  that  the  commonly  used 
linear  models  of  PEMFC  and  SOFC  are  not  controllable.  In  the  second 
part  of  the  paper  we  introduce  some  algebraic  constraints  on  these 
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models  that  follow  from  system  dynamic  equations  and  from  the 
steady  state  analysis. 

The  linear  mathematical  model,  for  the  PEMFC  dynamics  of  three 
fundamental  fuel  cells  dynamic  variables:  pressures  of  hydrogen, 
oxygen,  and  water  vapor,  was  derived  in  2004  [4].  The  model  was 
obtained  by  keeping  the  same  state  equation  and  slightly  modify¬ 
ing  the  output  equation  of  the  mathematical  model  derived  for  the 
SOFC  dynamics  in  2000  [5].  These  linear  mathematical  models  for 
PEMFC  and  SOFC  have  been  used  in  many  papers  including  some 
published  a  year  ago,  see  for  example  [6-13]  for  PEMFC  related 
problems  and  [14-16]  for  SOFC  related  problems. 

The  system  state  space  model  given  in  [4,5]  was  defined  by 
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with  the  state  space  variables  representing 

x(t)=[xi(t)  x2(t)  x3(t)]r  =  [pn2(t)  Po2(t)  PH2o(t)]r  (2) 

The  output  equation  represents  the  measured  fuel  cell  volt¬ 
age  and  it  is  obtained  using  the  Nernst  formula  for  the  open-loop 
cell  voltage,  V0(t),  and  subtracting  losses  due  to  the  cell  activation, 
Vact(t ),  and  due  to  the  stack  (fuel  cell)  resistance,  Vohm{t) 

ypEM(t)  =  VPEM(t)  =  V0(t)  -  VW(0  -  V0hm(t ) 

.N(£o+|ln{~!!}) 

-B  ln(C7(t))  -  RintI(t) 

The  system  inputs  are  molar  flow  rates  of  hydrogen  and  oxy¬ 
gen,  that  is,  qu2  (t )  and  qo2{t)  that  can  be  regulated  (controlled).  The 
stack  current  J(t)  plays  a  role  of  a  disturbance.  Note  that  C/(t)  must 
be  greater  than  1 ,  otherwise  the  activation  voltage  will  be  negative, 
and  hence  it  will  increase  the  open-loop  voltage  (instead  of  reduc¬ 
ing  it).  All  other  coefficients  are  assumed  to  be  constant.  The  values 
of  the  constant  coefficients  defined  in  the  model  equations  can  be 
found  in  [4]. 

The  SOFC  fuel  cell  model  of  [5]  has  exactly  the  same  state  Eqs. 
(1)  and  (2),  but  different  output  equation  in  which  the  activation 
voltage  is  not  present  in  the  expression  for  the  cell  output  voltage, 
that  is 

yso(t)  =  Vso(t)  =  JV  ^Eo  +  ^ln  j -  Rint/(t)  (4) 

Of  course,  in  the  SOFC  mathematical  model  (1),  (2)  and  (4)  the 
parameters  take  different  values  (except  for  the  universal  gas  con¬ 
stant  R  and  the  Faraday  constant  F).  The  values  of  the  constant 
parameters  for  the  SOFC  model  defined  by  (1 )  and  (4)  can  be  found 
in  [5], 

2.  Controllability  of  linear  PEMFC  and  SOFC  models 

The  importance  of  controllability  in  the  design  of  linear  con¬ 
trollers  for  PEM  fuel  cells  was  nicely  demonstrated  in  [17],  where 
even  for  originally  controllable  operating  points  of  a  linearized  sys¬ 
tem  some  design  techniques  provide  high  controllability  measures 
(requiring  less  control  efforts  and  more  efficient  control)  than  the 
other  also  controllable  operation  points.  Controllability  analysis  of 
liquid  water  in  a  fuel  cell  has  been  considered  in  a  very  recent  paper 
[18],  where  it  has  been  concluded  that  liquid  water  controllability 
is  needed  to  prevent  the  fuel  cell  flooding.  In  this  brief  note  we 
will  show  that  the  models  of  [4,5]  are  uncontrollable  (zero  control¬ 
lability  measure),  meaning  that  no  control  efforts  exist  to  satisfy 
general  goals  of  transferring  state  variable  from  a  given  initial  state 
to  a  desired  final  state  in  a  finite  time  interval  [3]. 

The  state  space  model  (1 )  can  be  represented  in  the  state  space 
form  as 
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where  u(t)  =  [  9^(0  9 o2 (0  1  are  the  control  inputs  and  d(t)  =  /(t) 
denotes  the  system  disturbance.  Using  the  standard  controllability 
test  [3],  we  can  form  the  controllability  matrix  for  the  state  space 
system  defined  in  (5)  given  by 

C(A,  B)  =  [  B  AB  A2B  ] 
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It  is  obvious  that  the  rank  of  the  controllability  matrix  C(A,  B )  is 
equal  to  2,  that  is 

rank{C(A,  B)}  =  2  <  3  =  n  (7) 

This  indicates  that  in  the  third-order  dimensional  linear  system 
considered  only  two  state  variables  are  controllable,  and  the  third 
one  is  uncontrollable.  Examining  the  state  space  equations,  it  can 
be  observed  that  the  equation  for  the  water  vapor  is  not  affected 
by  the  control  input  signal  and  hence,  the  water  vapor  pressure  is 
the  uncontrollable  variable  in  this  system. 

The  system  controllability  in  this  particular  model  is  needed  for 
several  reasons.  First  of  all,  it  is  well  known  that  being  uncontrol¬ 
lable,  the  state  variable  x3(t)  will  not  appear  in  the  system  transfer 
function  [3],  which  in  this  case  means  that  the  system  transfer 
function  is  of  order  two,  corresponding  to  the  controllable  state 
variables  Xi(t)  and  x2(t).  Flence,  every  frequency  domain  analysis 
that  involves  model  (1)  will  be  superficial  since  it  will  not  involve 
the  state  variable  x3(t).  Secondly,  it  is  known  from  [3]  that  state 
feedback  can  be  used  to  stabilize  unstable  systems,  but  it  can  not 
make  uncontrollable  systems  controllable  so  that  the  variable  x3(t) 
by  no  means  can  be  affected  by  control  input  signals,  and  it  will 
remain  affected  only  by  the  disturbance  signal  /(t)  that  changes 
randomly  as  /(t)  =  Vyc(t)/FL,  where  the  load  RL  changes  randomly 
in  time  as  a  piecewise  constant.  Hence,  changes  in  the  dynam¬ 
ics  of  the  state  variable  x3(t)  will  be  fully  determined  only  by  its 
time  constant  and  the  fuel  cell  disturbance  (current).  It  should  be 
emphasized  that  according  to  the  numerical  data  from  [4]  and  [5] 
the  time  constant  for  x3(t)  is  much  larger  than  for  the  remaining 
two  state  variables  (rH2o  =  18.418  s,  r02  =  6.64s,  rH2  =  3.37  s 
for  PEMFC  [4],  and  rH2o  =  78.3  s,  Tq2  =  2.91  s,  rH2  =  26.1  s  for 
SOFC  [5])  which  means  that  x3(t)  takes  much  longer  time  to  reach 
its  steady  state  value  (when  it  will  be  dictated  only  by  the  steady 
state  value  of  the  current  x^s  =  2/Cr/ss//CH2o)  than  the  remaining  two 
state  variables.  Moreover,  the  magnitude  of  the  state  variable  x3(t) 
is  much  smaller  than  the  magnitude  of  the  state  variables  X\  (t)  and 
x2(  t),  and  since  it  appears  in  the  dominator  of  the  cell  output  voltage 
formula  (4),  it  will  have  a  more  dominant,  more  lasting,  and  more 
unpredictable  impact  on  the  cell  output  voltage.  Thirdly,  control¬ 
ling  water  in  a  fuel  cell  is  fundamentally  important  [18],  since  it  can 
cause  cell  flooding,  degrade  the  cell  polarization  curve,  and  eventu¬ 
ally  damage  the  cell  membrane  [19,20]  (note  that  the  water  vapor 
mass  mH2o  is  proportional  to  the  water  vapor  pressure  pH2o  [19]). 

3.  System  analysis  constraints  of  PEMFC  and  SOFC  models 

In  this  section,  we  derive  some  algebraic  constraints  that  follow 
from  the  model  differential  equations.  These  constraints  were  not 
imposed  in  the  papers  [4,5]  that  derived  the  considered  models  nor 
in  any  other  follow-up  paper  that  have  used  these  models  alone 
or  in  hybrid  configurations  with  other  electric  energy  generating 
devices.  The  constraints  are  imposed  at  steady  state,  for  the  initial 
conditions,  and  for  all  time  instants. 
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3.1.  Steady  state  constraints 


The  steady  state  value  of  the  water  vapor  (obtained  by  set¬ 
ting  the  derivative  to  zero  and  solving  the  corresponding  algebraic 
equation)  is  given  by 
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The  steady  state  values  of  the  hydrogen  and  oxygen  pressures 
are  functions  of  the  steady  state  values  of  the  molar  inlet  rates  of 
hydrogen  and  oxygen  and  the  steady  state  fuel  cell  current.  They 
are  respectively  given  by 
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The  corresponding  steady  state  output  voltage  value  (for  SOFC) 
from  (4)  is  given  by 
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Since  the  hydrogen  and  oxygen  pressures  are  positive  quanti¬ 
ties,  the  following  conditions  on  the  inlet  hydrogen  and  oxygen 
flow  rates  must  be  satisfied  at  the  steady  state 
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Another  set  of  simulation  constraints  comes  from  the  Nernst 
open-loop  voltage  that  is  a  positive  quantity.  Since  the  “In”  opera¬ 
tion  is  present  in  that  formula  we  must  have 


The  analytical  expression  for  the  hydrogen  and  oxygen  pres¬ 
sures  and  corresponding  initial  condition  constraints  are  given  by 
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It  is  interesting  to  observe  the  initial  condition  constraints  on 
xi(0)  and  x2(0)  imposed  in  (15)  and  (16).  Even  though,  the  coeffi¬ 
cients  that  multiply  the  stack  current  in  (15)  and  (16)  are  small  (in 
[4]  they  are  0.02  and  0.01,  respectively  for  hydrogen  and  oxygen) 
these  initial  condition  constraints  must  be  satisfied  in  simulation 
studies  whenever  these  models  are  used,  otherwise  in  the  initial 
time  interval  the  negative  values  could  be  obtained  for  the  hydro¬ 
gen  and/or  oxygen  pressures  (using  MATLAB/Simulink  or  any  other 
simulation  software  package).  There  is  no  similar  initial  condition 
constraint  on  the  water  vapor,  except  for  the  obvious  one  that 
comes  from  the  Nernst  formula  x3(0)  ^  0  (also,  the  same  formula 
require  in  general  x3(t)  ^  0,  V  t).  The  water  vapor  values  in  time  as 
a  function  of  the  fuel  cell  current  are  obtained  by  solving  the  last 
equation  in  (1 ),  leading  to 
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clearly  indicating  that  the  expression  for  the  water  pressure  is 
not  a  function  of  the  inlet  molar  flow  rates  of  the  hydrogen  and 
oxygen. 

Since  the  cell  current  depends  on  the  load  RL  that  changes  in 
time  independently  outside  of  the  fuel  cell  (as  a  piecewise  constant 
function),  with  Vjc(t)  =  i?J(t),  the  cell  current  is  in  general  variable 
and  the  above  stated  constraints,  (11)  and  (13)-(16)  must  hold  for 
the  worse  case  scenario  (when  J(t)  takes  its  maximal  values,  that  is, 
for  /max  =  max  {/(t)} ,  V  t). 
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Hence,  the  steady  state  hydrogen  and  oxygen  pressures  must 
be  such  to  overcome  not  just  the  constraint  given  in  (1 1 ),  but  also 
the  stronger  constraint  given  in  (13),  which  is  dictated  by  the  cell 
(stack)  steady  state  current. 


4.  Conclusions 

The  lack  of  controllability  of  considered  linear  models  of  PEM 
and  SO  fuel  cells  might  mean  that  the  results  presented  in  journal 
and  conference  papers  using  these  uncontrollable  models  might  be 
valid  only  for  the  chosen  set  of  data,  and  that  the  conclusions  drawn 
in  those  papers  are  not  general  (valid  for  all  possible  inputs  and  all 
possible  values  of  the  state  space  variables).  We  have  also  intro¬ 
duced  some  algebraic  constraints  that  must  be  satisfied  at  steady 
state,  initial  time,  and  at  all  times  to  be  able  to  run  simulations  (in 
MATLAB/Simulink  or  any  other  computer  software  package)  of  the 
considered  linear  PEM  and  SO  fuel  cell  models. 


3.2.  Initial  condition  and  time  constraints 
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The  pressure  positivity  requirements  will  also  require  imposing 
constraints  on  the  hydrogen  and  oxygen  pressure  initial  conditions. 
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